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Design and synthesis of two-dimensional (2D) materials with robust ferromagnetism and biocompatibility 
is highly desirable due to their potential applications in spintronics and biodevices. However, the hotly 
pursued 2D sheets including pristine graphene, monolayer BN, and layered transition metal 
dichalcogenides are nonmagnetic or weakly magnetic. Using biomimetic particle swarm optimization 
(PSO) technique combined with ab initio calculations we predict the existence of a 2D structure, a 
monolayer of rocksalt- structured CrN (100) surface, which is both ferromagnetic and biocompatible. Its 
dynamic, thermal and magnetic stabilities are confirmed by carrying out a variety of state-of-the-art 
theoretical calculations. Analyses of its band structure and density of states reveal that this material is 
half-metallic, and the origin of the ferromagnetism is due to p-d exchange interaction between the Cr and N 
atoms. We demonstrate that the displayed ferromagnetism is robust against thermal and mechanical 
perturbations. The corresponding Curie temperature is about 675 K which is higher than that of most 
previously studied 2D monolayers. 



Two dimensional (2D) materials have received intensive attention^"^ owing to the birth of graphene^. 
Materials in low dimensions can have electronic and magnetic properties that are fundamentally different 
from those of bulk phases. For instance, cubic boron nitride (BN) is an insulator whereas BN nanofilms can 
be intrinsically metallic^ Bulk M0S2 is an indirect band gap semiconductor while monolayer M0S2 has a direct 
bandgap of 1.90 eV, making it a promising material for nanoelectronics and valleytronics^'^. In this work, we 
explore the possibility of achieving highly spin-polarized ferromagnetism in chromium nitride (CrN) by altering 
its dimensionality. 

Transition metal nitrides (TMNs) constitute a large family of materials with rich fundamental properties and 
vast technological applications. They have traditionally been used at extreme conditions of temperature and 
pressure due to their superhardness, high melting points, and outstanding corrosion resistance. In this family, 
CrN, a prototypical TMN, is particularly interesting not only because it combines all the favorable properties of 
TMNs but also due to its high thermal stability, good electrical conductivity^, and especially its biocompatibility 
giving rise to its applications on medical implants and devices, such as artificial joints and cutting tools^. CrN has 
challenged both theoretical and experimental investigations^" for its rare "itinerant antiferromagnetism"^\ 
Currently, two typical phases have been identified: the rocksalt- structured phase with Fm-3m symmetry and 
paramagnetic (PM) order, and the orthogonal phase with Pnma symmetry and antiferromagnetic (AFM) 
Qj-jgj-8,11 Driven by magnetic stress induced lattice distortion, the former phase can transform into the latter 
at a very low external pressure^^'^^. However, neither the rocksalt-structured nor the orthogonal phase exhibits net 
magnetization despite remarkable local magnetic moments observed. 

Previous studies have demonstrated that magnetism of nanomaterials differs from that of bulk materials in a 
fundamental way because both magnetic coupling and magnetic moments are strongly dependent on coordina- 
tion number and dimensionality- related quantum confinements. For example, in nanoscale, the AFM coupling 
between Cr atoms can be switched to ferromagnetic (FM) by means of introducing nitrogen^^ or oxygen atoms^^. 
The magnetic coupling between the Mn/Cu dopants in dilute magnetic semiconductors (DMSs) is highly 
dependent on dimensionality, namely, exhibiting FM coupling configuration in ID GaN nanowires^°"^^ while 
having AFM coupling in the 2D thin films^^, which implies the possibility of achieving ferromagnetism in CrN 
system by altering its dimensionality. What's more, even in the AFM phase of bulk CrN there exists ferromagnetic 
layers^' ^\ These inspired us to carry out the following study to seek a FM CrN system as it may have applications in 
spintronics as well as biotechnology. 
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Figure 1 | CrN in 2D form. The 2D structures of CrN predicted by 
CALYPSO with respect to the number of formula units Z per supercell. 
Silver and blue balls represent Cr and N ions, respectively. For clarity all 
structures are shown in 2 X 2 supercells (green lines). The red lines denote 
the lattice basis vectors of the unit cell. 

Results 

Structure search. We first revisit CrN in bulk phase to test the 
methods that we used. Full relaxations of the lattices and atomic 
positions are performed for the rocksalt phase of CrN with both 
the FM and AFM coupling configurations of the Cr atoms, 
respectively. The FM configuration is found to be energetically less 
favorable than the AFM one by an energy difference of 210 meV, 
which is in good agreement with previous results^ as well as the fact 
that bulk CrN is an AFM material at low temperature. We then begin 
the study of CrN in two-dimensional form. A thorough exploration 
on the potential energy surface (PES) of 2D CrN is needed to identify 
the thermodynamically stable structure. We perform unbiased 
structural searches by means of the CALYPSO (Crystal structure 
AnaLYsisby Particle Swarm Optimization)^^ technique based on 
the biomimetic Particle Swarm Optimization (PSO) algorithm, 
which has successfully predicted a series of 2D crystals^^"^^ with a 
considerably high efficiency. Without any priori structural 
information, our blind structure search converges to a planar one- 
atomic- thick structure with an optimized lattice constant of 2.82 A 
and Cr-N bond length of 1.99 A, irrespective to the choice of 
supercell size Z (Z represent the number of formula units in a 
supercell,) as shown in Figure 1 (See the Supplementary 
Information for details on the structure search). 

Preferred Magnetic Configuration. The predicted 2D CrN 
structure turns out to be surprisingly simple, which can be viewed 
as a monolayer of the rocksalt-structured CrN(lOO) surface. To 
search for the magnetic ground state of this monoatomic thick 
structure, we calculate the spin polarization energy^^, namely, the 
energy difference between the nonmagnetic and magnetic 
configurations. Our calculations yield a large energy gain of 
1.6 eV/unit cell upon spin polarization, indicating that the 
magnetic configuration is highly preferred and the monolayer CrN 
might be a good magnetic material. In the spin polarized 
configuration, the bond length of Cr-N is 1.99 A, which is 4.7% 
larger than that of the nonmagnetic configuration. Each unit cell 
carries a net magnetic moment of 3.00 |iB, which arises mainly 
from the Cr 3d orbitals (3.19 jis) and the neighboring N atom is 
polarized antiferromagnetically with a magnetic moment of —0.37 
jiB- This is because the Cr ions in this nanosheet have valance state of 
+ III based on the chemical stoichiometry and the electronic 
configuration of Cr^^ is 3cP. We then consider a low spin 
configuration by constraining the total magnetic moment of a unit 
cell to be 1 ^b- The calculated total energy of the low spin 
configuration is 0.95 eV higher than that of the high-spin state. In 
addition, we find the lattice constants undergo a expansion from 
2.70 A X 2.70 A (the NM state) to 2.73 A X 2.75 A (low-spin state) 
to 2.82 A X 2.82 A (high-spin state), and the Cr-N bond length changes 
from 1.90 A to 1.94 A to 1.99 A, accordingly. This indicates that spin 
polarization can effectively release the strain and stabilize the system. 
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Figure 2 | Magnetic coupling configurations of the monolayer CrN. Top 

and side views of the optimized geometry and spin charge density 
isosurface (isovalue: 0.2 e/A^) of (a) the FM and (b) AFM configurations of 
the monolayer CrN. Red and yellow isosurfaces correspond to spin up and 
spin down charge density, respectively. Green squares denote the 2X2 
supercell. 

To further determine the most preferable magnetic coupling 
between the Cr atoms in this predicted CrN nanosheet, we perform 
spin polarized total energy calculations for the two typical collinear 
magnetic configurations, namely, the FM and AFM states, within a 2 
X 2 supercell. The FM state is found to be 632 meV lower in energy 
than the AFM one, indicating that the monolayer CrN is FM. Such a 
magnitude of energy difference is observable in experiments. Figure 2 
shows the spin charge density distributions (p|-p|) for both the FM 
and AFM coupling configurations. We note that large spin charge 
accumulates around the Cr^^ sites in both the two cases, resulting in 
large local magnetic moments. However, the spin charge distribution 
of the N 2p electrons appears differently: In the FM ground state, the 
N 2p electrons are strongly spin polarized antiferromagnetically with 
the Cr ions, while in the AFM configuration the net spin polarization 
in the ligands is negligible, resembling that in its bulk counterpart^^. 

Dynamic Stability. To ensure that the CrN nanosheet can be 
practically utilized as a FM material, we first examine its 
dynamical stability by calculating its phonon spectra. The Phonon 
dispersion and frequency density of states (PhDOS) are displayed in 
Figure 3(a). All vibrational modes are found to be real in the whole 
Brillouin Zone, indicating this structure is dynamically stable and 
can exist as a freestanding 2D crystal. 

Thermal Stability. We then examine the thermal stability of the 
monolayer CrN. We construct a large supercell of 6 X 6 
containing 36 formula units to minimize the constraint induced by 
periodicity. Ab initio molecular dynamics (MD) simulations^^"^^ are 
performed at 300 K with a time step of 1 fs using Nose heat bath 
scheme^^ to investigate whether the FM monolayer is thermally 
stable at room temperature. The atomic configuration of the 
monolayer remains nearly intact after heated for more than 6 
picoseconds. The variations of the total potential energy and total 
magnetic moment with respect to simulation time are plotted in 
Figure 3(b), showing that the monolayer is thermally stable as both 
of them remain almost invariant during the simulation. We further 
gradually increase heating temperature and find the structure can 
survive up to 600 K. MD simulations are also carried out for a 4^2 X 
4^2 supercell containing 32 formula units. We find, again, neither 
structural reconstruction nor magnetic phase transition occurs when 
the structure is heated at room temperature (see Figure S2 in 
Supplementary Information for details). These results demonstrate 
that the FM monolayer CrN, once synthesized, is stable and can stay 
robust at room temperature. 
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Figure 3 | Structural and magnetic stability of the monolayer CrN. (a) Phonon dispersion and frequency DOS (in arbitrary unit) of the monolayer 
CrN.The high symmetry q point path in the Brillouin Zone is r(0,0) X ( 1/2, 0) ^ M (1/2, 1/2) Y (0, 0). (b) Variations of the total potential energy 
and total magnetic moment with respect to simulation time. ( c) Strain energy of the FM and AFM configurations under the biaxial strains. The inset shows 
the 2X2 supercell used to calculate the strain energy, (d) Variations of the average magnetic moment per unit cell and magnetic susceptibility {%) with 
respect to temperature. 



To check if any possible structure reconstruction or buckling pat- 
tern^^ could happen, we create a single-layer PbO-type structured 
CrN which is ubiquitous in many layered crystals^^. We find, upon 
full relaxation, the PbO-type structure of CrN spontaneously con- 
verges to the planar monolayer sheet, indicating the outstanding 
stability of this CrN structure. 

Magnetic Stability. We further study the stability of the FM order of 
the monolayer CrN against the in-plane strains. CrN thin films, in 
reality, usually grow on certain substrates like MgO (100) surface^'^^. 
The synthesized thin films might endure some external strain, either 
compressive or tensile stain, due to the possible lattice mismatch. We 
calculate the strain energy for both the FM and AFM configurations 
under the isotropic biaxial in -plane strain. The results are plotted in 
Figure 3(c). Although the exchange energy decreases with the 
compressive biaxial strain and increases with the tensile strain, the 
FM state is found to always lie lower in energy than the AFM 
configuration in a wide range of the biaxial strains from —6% to 
8%. This indicates that the FM order in the monolayer CrN can 
remain virtually undisturbed from the in-plane strains. A 
noteworthy trend in the strain energy curves is that the FM state 



becomes less competitive upon increasing the compressive biaxial 
strain, which is reminiscent of previous prediction on the bond 
length contraction induced AFM ordering in the Mn-doped GaN 
thin films^^. 

Estimation of Curie temperature. Having confirmed the structural 
and magnetic stability, we then quantitatively study the stability of 
the FM ordering against thermal perturbation by estimating its Curie 
temperature. To this end we carry out Metropolis Monte Carlo (MC) 
simulations based on 2D Ising modeP. According to the classical 
Heisenberg model, the Hamiltonian of a 2D lattice can be expressed 
as H = — ^^^Jrhirhj, where m/ corresponds to the magnetic 

moments (in ^b) at site / (each unit cell is treated as one site, and 
for each site four spin states are considered, namely, two high spin 
sates, 8=3/2 and S=— 3/2, and two low spin states, 8 = 1/2 and 
8 = — 1/2, and / refers to the Heisenberg direct exchange parameter, 
which can be obtained by / = {l/S)EeJ2m^, where £ex represents the 
exchange energy of the AFM and FM states in the 2X2 supercell. 
The factors of 1/8 and 1/2 correspond to the eight exchange pairs in 
the supercell and the double count of each pair in the summation. 
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Figure 4 | Electronic properties. Spin-polarized electronic band structure (Red: spin up channel; Blue: spin down channel) and PDOS of the monolayer 
CrN. The Fermi level is shifted to 0 eV. 



respectively. Substituting E^^ = 632 meV and m = 3 into the 
formula, we get / = 4.39 meV. For simplicity, only the nearest 
neighbor exchange interaction is taken into account in the 
Hamiltonian. In order to reduce the periodic constraints, a 50 X 
50 supercell is used to mimic the monolayer sheet. The variations 
of average magnetic moment per unit cell and magnetic susceptibility 
with respect to temperature are plotted in Figure 3(d). We note that 
the magnetic moment remains in the high spin state in the low- 
temperature range and then drastically drops to near zero and the 
magnetic susceptibility appears a sharp peak at a critical temperature, 
showing a thermal-induced FM-to-PM phase transition. The Curie 
temperature is estimated to be around 675 K, which is higher than 
that of most previously studied 2D systems^^'^^'^^"^^. This implies that 
ferromagnetism in the monolayer CrN is robust and can be detected 
at room temperature. 

Electronic Properties. Bulk CrN is known as a "borderline Mott 
insulator" to investigate the electronic properties of the 
monolayer sheet, we calculate the spin-polarized electronic band 
structure. The results are plotted in Figure 4. We note although the 
band dispersion relationships in both spin channels are similar, it 
exhibits metallic feature in the spin-up channel characterized by the 
Fermi level crossing the partially occupied valance bands in the 
Brillouin zone and semiconducting feature in spin-down channel, 
indicating that the monolayer CrN is half-metallic. We further 
calculate the total density of states (DOS) and orbital resolved 
partial DOS at the Cr and N sites to explore the origin of the half- 
metallicity and ferromagnetism. There is a remarkable DOS peak at 
the Fermi level which is only contributed by the spin up channel. The 
partial DOS shows an explicit spin splitting in both the Cr 3d and 
N 2p states. The overlap between the Cr 3d and N 2p orbitals leads to 
the significant DOS near the Fermi level in spin-up channel, thus 
results in the half-metallicity as well as the ferromagnetism. 
Observably, the N 2p orbitals play key role in mediating the 
magnetic coupling between the Cr atoms. 

Bilayer Stacking: structural and magnetic configurations. So far, 
we have established the viability of realizing robust ferromagnetism 
in the monolayer CrN. It would be intriguing to see how the 
structural and magnetic configurations change when the two CrN 



monolayers are stacked to form bilayer CrN. To address this concern 
we build two typical bilayer models, namely, AA and AB stacking, as 
shown in Figure 5. Our calculated results show that the AB stacking is 
energetically more favorable than the AA stacking. The binding 
energy between the two monolayers is 1.428 eV and 0.641 eV, and 
the optimized interlayer distance is 2.09 A and 3.18 A for the AB and 
AA stacking patterns, respectively, showing different bonding 
characters, namely, weak interaction in the AA stacking and a 
strong bonding with a hybrid ionic and covalent feature in the AB 
stacking. This is different from the situations of graphene and BN 
sheets which have weak interactions between layers, but is similar to 
the case of Fe-phthalocyanine sheets, where the interaction between 
the two layers is strong. However, a single Fe-phthalocyanine sheet 
can be synthesized on an insulating substrate^^, where the Fe- 
phthalocyanine sheet only weakly interacts with the substrate^^. 
Therefore, we can expect that our predicted CrN sheet could be 
grown on some insulating substrates. Since ferromagnetism of the 
monolayer CrN is energetically robust against antiferromagnetism, 
for simplicity we first assume that the intra-layer coupling is FM like 
a freestanding monolayer, and for each stacking pattern we consider 
two inter-layer couplings, namely, FM and AFM. We relax the lattice 
and atomic positions with inclusion of Grimme's potentiaP° to 
compensate for missing van der Waals interactions in standard 




fex = -96 meV ^ex = 209 meV 

(a) (b) 

Figure 5 | Stacking pattern of Bilayer CrN. (a) AA stacking and (b) AB 
stacking. The exchange energies (define as energy difference between the 
AFM and FM configurations) of possible magnetic configurations are 
presented. Here FM and AFM refer to inter-layer magnetic coupling. 
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DFT calculations. Total energy calculations show that the AB 
stacking FM coupling configuration has the lowest energy among 
all the considered bilayer models. We then enlarge the AB stacking 
model to a 2 X 2 supercell to examine the intra-layer magnetic 
coupling. Three possible magnetic coupling configurations: FM 
(both intra- and inter layer couplings are FM), partial AFM (intra- 
layer FM and inter-layer AFM), and AFM (both intra- and inter- 
layer couplings are AFM) are considered. The calculated total energy 
relative to the FM state for the partial AFM and AFM state are 
248 meV and 209 meV, respectively. This confirms that the 
bilayer CrN is also FM. Therefore, we can conclude that the CrN 
monolayer, either in freestanding form or stacked in energetically 
favorable patterns, can exhibit robust ferromagnetism. 

Discussion 

To better understand the ferromagnetism of the monolayer CrN, it is 
worth recalling the following facts. (1) Chromium, due to its unique 
electronic configuration of [Ar]3(i^4s^ and the small energy gap 
between the 3d and 4s orbitals, can exist in different oxidation states, 
namely, 0, +2, +3, and +6, which can be used to tune the structural 
and magnetic properties of Cr- containing materials by changing its 
bonding feature. (2) In a Cr2 dimer, Cr atom binds strongly with the 
other Cr atom and the resulting sextuple bonding yields a large 
binding energy (1.44 eV) and a very short bond length (1.68 A). 
Small clusters of Cr containing five atoms or fewer are all AFM^\ 
In the bulk phase, Cr is also AFM^^. However, the structures and 
magnetic properties of these small clusters can be drastically modi- 
fied by doping a nitrogen atom^^, where the nearest-neighbor Cr 
atoms are coupled antiferromagnetically to the N atom and hence 
ferromagnetically with each other. Similarly, in the ground state of 
the monolayer CrN, the two nearest Cr atoms are too far apart to 
allow significant direct overlap of the unpaired electron wave func- 
tions, the coupling involves the N 2p orbitals polarized antiferromag- 
netically to the 3d orbitals of Cr atoms. This lowers the energy of the 
system and makes the Cr atoms couple ferromagnetically to each 
other. Thus, the mechanism responsible for the FM coupling in the 
monolayer CrN can be classified as the p-d exchange mechanism^^. 
(3) To achieve highly spin polarized ferromagnetism, usually it is 
necessary to introduce transition metal atoms. The observed mag- 
netism originated from 2p electrons in the currently studied systems 
like porous graphene or hydrogenated^^ or fluorinated graphene is 
weak. Among 3d transition metals, Cr and Mn have been widely used 
as dopants in semiconductors such as GaN or GaAs or ZnO to 
introduce magnetism due to their unique electronic configurations. 
However, it is challenging to regularly dock Cr or Mn atoms into the 
bulk phases and even on the 2D structures because of the clustering 
of TM atoms resulting from strong d-d orbital interactions. Recent 
experimental advances have demonstrated that 2D phthalocyanine- 
based organometallic (labeled as TM-Pc) sheets^°'^^ and tetragonal 
transition-metal-based 7,7,8,8-tetracyanoquinodimethane molecule 
coordination sheets (labeled as TM@TCNQ)^^ have been success- 
fully synthesized with the TM atoms regularly and separately dis- 
tributed, providing the possibility of synthesizing and the 2D CrN 
sheet experimentally and detecting the ferromagnetism. (4) Different 
from Mn-containing 2D sheets where Mn^^ and Mn^^are neuro- 
toxic^^, as they can enter the brain by crossing the blood— brain 
barrier and cause permanent neurological disorder, the explored 
2D CrN sheet is biocompatible, where trivalent chromium (Cr^^) 
in CrN is biologically active^^, widely found in food and supplements, 
meanwhile nitrogen can also improve biocompatibility significantly^^. 

Summary 

Motivated by design of novel 2D magnetic materials and inspired by 
the chemical control of magnetic coupling in molecules and clus- 
ters^^'^^, a comprehensive theoretical study on geometry, stability, 
electronic properties and ferromagnetism of the monolayer CrN 



sheet is performed. We find CrN, a well-known AFM or PM material 
in bulk phases, can become spontaneously FM when sUced into one- 
atom-thick nanosheet. The switching between AFM or PM to FM 
couplings is brought by tailoring the dimensionality and the origin of 
ferromagnetism is ascribed to p-d exchange mechanism. The well- 
defined atomic configuration and intrinsic ferromagnetism of the 
monolayer CrN suggest it is superior to other 2D magnetic materials: 
(1) Unlike other 2D Cr-containing sheets such as Cr-Pc and TM@ 
TCNQ monolayers, which are AFM, the CrN nanosheet predicted 
here is both stable and FM, and the ferromagnetism can be sustained 
against mechanical and thermal perturbations; (2) Unlike other TM- 
based 2D materials which needs external stimulus like chemical 
functionalization^^'^^'^^, carrier injection^^'*^ or external strain^^ to 
switch on or enhance the FM order, the 2D CrN sheet possesses 
intrinsic ferromagnetism with a very high Curie temperature of 
675 K, which is much higher than that of previously reported 
150 K and 140 K for Mn-Pc sheet^° and Mn02 monolayer^^, respect- 
ively; (3) Unlike conventional Cr or Mn doped DMSs where semi- 
conductors are used as matrix to host the magnetic dopants, and 
magnetism, in most cases, comes from the clustering of TM atoms, 
in the monolayer CrN, Cr ions are individually held at the specific 
sites in a 2D crystal structure, and couple instrically ferromagneti- 
cally with each other. The Curie temperature is much higher than the 
highest Tc value achieved experimentally in Mn-doped GaAs. In 
addition, CrN thin film has been widely used as a coating material 
owing to its outstanding corrosion-resistance^". Therefore, the dis- 
played high stability, robust ferromagnetism and the intriguing bio- 
compatibility make this monolayer CrN sheet a very promising 
candidate for applications in spintronics, biodevices and functional 
coating techniques. We hope that the present study will stimulate 
further experimental effort in this direction. 

Methods 

Our first principles calculations are based on spin-polarized Density Functional 
Theory (DFT). All-electron Projector Augmented Wave (PAW) method imple- 
mented in Vienna Ab initio Simulation Package (VASP)^'' is used to investigate the 
ground state structural, electronic and magnetic properties of CrN in the bulk and 
monolayer phases. For N and Cr, 2s^2p^ and 3c/^4s^ electrons are treated as valence 
electrons, respectively. Considering the localization character of Cr 3c/ electrons, 
GGA+ U formalism is adopted to describe the strong on-site Coulomb interaction. 
We choose the simplified (rotationally invariant) approach due to Dudarev et aV^, 
which involves only one parameter, namely, the effective U{U-J, where 17 is the onsite 
Coulomb repulsion and / is the onsite Hund exchange coupling). We set the effective 
(7 value to be 3 eV, which has been widely used in both bulk^^'^^'^^ *^^ and low- 
dimensionaP"''*'^^'*^* systems that contain Cr or other 3d transition metal atoms and 
provided reasonable results. 2D structure searches are carried out using the 
CALYPSO methodology^"^"^^ in combination with VASP for local optimizations and 
total energy calculations. Phonon properties are calculated using Density Functional 
Perturbation Theory (DFPT)''' implemented in the Quantum ESPRESSO suite''^ Ab 
Initio Molecular Dynamics (AIMD) simulations are performed to examine the 
thermal stability. Detailed computational procedures are provided in the 
Supplementary Information. 
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